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Q ■ Abstract 

O . We report on studies of multi-parton corrections from nonlocal operator ex- 

pansion. We discuss relations between eikonal-line matrix elements and parton 
distributions, and present an illustration for initial-state coUinear evolution. 
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1 Introduction 

^ _ hadronization, soft underlying scattering, multiple hard interactions. Models for these processes are 

Q ' necessary, for instance, for Monte Carlo generators to produce realistic event simulations. 

' The treatment of multiple parton interactions in QCD will require methods that go beyond the 

local operator expansion, and likely involve fully unintegrated parton correlation functions [1]. Besides 

, the relevance for event generators, this should also provides a natural framework for the investigation at 

(yr^ ' the LHC of possible new strong-interaction effects at very high energies, including parton saturation [2]. 
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This report is based on the analysis [3] of nonlocal operator expansion, investigating corrections 
^ ■ from graphs with multiple gluon exchange. The point of view in this study is to connect the treatment of 

^ . multi-gluon contributions with formulations in terms of standard partonic operators, and in this respect 

it can be seen as deriving from the approach of [4]. We present an illustration for the case of structure 
functions. This case is also treated in the analyses of [5]. More discussion may be found in [6]. The 
formulation discussed below trades parton distribution functions for moments of eikonal-line correlators. 
We expect this formulation to be useful also for the treatment of the associated final-state distributions. 
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I 2 From parton distribution functions to eikonal-line matrix elements 

The analysis [3] starts with the quark distribution function, defined as 

f,{x,^^) = dy-e'-P^y- {P\mQ{d)^+Q\y-)mv- MP)c (D 

where ifj is the quark field, Q is the gauge link, and the subscript c is the instruction to take connected 
graphs. The matrix element ([T]) can be rewritten as the real part of a forward scattering amplitude [3], in 
which we think of the operator Q'^ifj as creating an antiquark plus an eikonal line in the minus direction, 
starting at distance y^ from the position of the target. 

Next, supposing that x is small, we treat the evolution of the antiquark-eikonal system in a hamil- 
tonian framework (see [3] and references therein) which allows us to express the evolution operator in 
the high-energy approximation as an expansion in Wilson-line matrix elements. The leading term of this 
is ("dipole" term) 

H(z,&) = j[dP']{P'\j^TT{l-F\h + z/2)F{b-z/2)}\P) , (2) 



where F is the eikonal operator 

F{r) =Vexp!^-ig J dz-^+(0, z", r) t^j , 

z is the transverse separation between the eikonals in and b is the impact parameter. 

In this representation the quark distribution ([Hi is given by the coordinate-space convolution 



(3) 



xfq{x,fj,) = J dh dz u{fi,z) E{z,b) — UV 



(4) 



In [3] the expUcit result is given for the function u{fi, z) at one loop in dimensional regularization and 
for the counterterm —UV of MS renormalization. The MS result can also be recast in a physically more 
transparent form in terms of a cut-off on the z integration region, as long as the scale /i is sufficiently 
large compared to the inverse hadron radius: 
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e{z^l? > a^)E{z,b) , 
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(5) 



where a is a renormalization scheme dependent coefficient given in [3]. 

The Wilson-line matrix element E{z,b) receives contribution from both long distances and short 
distances. At small z it may be treated by a short distance expansion. At large z it should be parameter- 
ized consistently with bounds from unitarity and saturation [2] and determined from data. 

3 An algebraic relation for eikonal operators 

A general relation between fundamental and adjoint representation for E, valid for any distance z, is 
given in [3], based on the algebraic relation 
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(6) 



withF = Ffund.,f/ = ^adj.- 

From this one can obtain small- 2; relations connecting S to the gluon distribution. For instance, 
for the fundamental representation at small z this yields 



E{b, z) = z^ xG{x, fj,) (j){b), 



(V) 



where by xG we denote the gluon distribution (either the Xc-scale or weighted-average expressions 
in [3]), and </>(&) obeys 

j db (l){b) = 1. (8) 

The result for H in the fundamental representation corresponds directly to the one for the dipole cross 
section in the saturation model [2]. Results in the fundamental and adjoint cases are relevant to discuss 
quark saturation and gluon saturation. 



4 Power-suppressed contributions 

In the s-channel framework of [3] contributions to hard processes suppressed by powers of the hard scale 
are controlled by moments of H, 

p f dz 



M, 



T{l-p) J TTZ 



(9) 



analytically continued for p > 1. Models for the dipole scattering function including saturation are 
reviewed in [2]. In this case the moments Q are proportional to integrals over impact parameter of 
powers of the saturation scale. Higher moments are obtained from derivatives with respect to p, 



db [Qlib)]P , Mj,,k ^ (-1)' ^^Pfi ■ (10) 

As an illustration, we determine the Ca/x part of the coefficients of the first subleading power 
correction from the s-channel for transverse and longitudinal structure functions Ft, F^. Denoting the 
derivative by Fj = dFj/dln for j = T, L, and its leading-power contribution by Fj iead.^ one has 

Fj - Fj^lead. = bj,o -M2,0/(3' + M2,l/Q^ + ■■■ (H) 

Structure functions can be analyzed in the same way [3] as described in Sec. [2] for the quark 
distribution function. The main difference compared to the case of the quark distribution ([T]) is that the 
ultraviolet region of small z is now regulated by the physical scale rather than requiring, e.g., MS 
renormalization. Saturation is reobtained [3] within the dipole approximation [2]. By the analysis based 
on I©,© the saturation scale Qs{b) for a dipole in the fundamental representation is 

Ql{b)='^^xG{x,^^)m■ (12) 

To study the expansion in powers ofl/Q^ it is convenient to go to Mellin moment space by representing 
E via the Mellin transform 

E{z,b) = z' / —{z')-^E{u,h) , (13) 



< a < 1. Then the structure functions Ft,l have the representation 

f a+ioo 



r ra+ioo _ 

xFt,l = db — H(n,b) ^t,l{u) , 

J Ja—ioo '^^'^ 



(14) 



where ^T,Liu) can be read from [7] and are given by 



, o^ Nr , ^n^„T(3 — u)T(2 — u)T(l — u) , ^ ^ 
*T W = (4>^(Q^)" V(5/2 -\.)r(3/2Vu) + ") "^W' ('« 

^.(") = <-^>iJ^('3^)° r(5/2'-n)'r(3T2 + u, ' ^" " "'• <'« 

with r the Euler gamma function. The expansion m\/Q^ of ([T4l) is controlled by the singularity structure 
of the integrand in the u-plane [3,5,6]. Eqs. ([T5]).(fT6l) show that longitudinal has no pole at u = 0, 
so that the leading singularity is given by the u = pole in H, while the first subleading pole u = — 1 is 
absent in transverse <I>t due to the numerator factor (1 + u), so that the answer for the transverse case at 
next-to-leading level is determined by the singularity in H, with $ contributing to the coefficient of the 
residue. 

It can be verified that contributions to (fT4l) in the lowest p = 1 moments in Eq. (ITOl ) correctly 
reproduce the small-x gluon part of renormalization-group evolution, 

■ / 2\ «s /"^(iz [z^ + (1 — z)^] . fx ^\ , 
FxMad. = {ej ^ J ~ 2 ■'^Vz' V + ^^^^^ 

— i^a) 77" TT + quark term , (17) 
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using ((8]). (fT2l) and the gluon distribution G evaluated at the average [3] x ~ Xc, with the lowest x-moment 
of the gluon quark splitting function 

f dz Pgg{z) = [ dz + (1 - zf]/2 = 1/3 . (18) 
Jo Jo 

Beyond leading power, the first subleading corrections read 

Ft - FT,iead. = -i^l)^^ I db [Qlib)? + • ■ ■ , (19) 

+ {el)-^^ldb[Ql{b)]Hn[Q'/Q'M] + --- ■ (20) 

That is, the b coefficients in (fTTI ) are given by 

bT,o = -{el) Ca/{20tt^x) , 6t,i = , 
&L,o = -(e^)CA [14/225 + V'(l)/15]/(7r3x) , bL,i = (el) CA/{l5Tr^x) , (21) 

with ip the Euler psi function. 

Via process-dependent coefficients analogous to those in (ITTI ). the eikonal-operator moments Q 
will also control power-like contributions to the associated jet cross sections due to multi-parton inter- 
actions in the initial state. At present these processes are modeled by Monte Carlo, which point to their 
quantitative significance for the proper simulation of hard events at the LHC. The above discussion sug- 
gests the potential usefulness in this context of analyzing jet and structure function data by trading parton 
distribution functions for s-channel correlators defined according to the method of Sec. |2] 
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